Mismatch Repair Error Implies Chargaff's Second Parity Rule 

Bo DengS 

Abstract: Chargaff's second parity rule (PR2) holds empirically for most types of DNA that 
along single strands of DNA the base contents are equal for complimentary bases, A = T, G = 
C. A Markov chain model is constructed to track the evolution of any single base position 
on a given single strand of DNA whose organism is equipped with the process of mismatch 
repair. Under the key assumptions that the mismatch error rates primarily depend the com- 
O ■ plementarity and the steric effect of the nucleotides and that the mismatch repairing process 

^ ■ itself makes strand recognition errors, the model shows that the steady state probabilities for 

p the base position to take on one of the 4 nucleotide bases are equal for complimentary bases. 

\ As a result, Chargaff's second parity rule is the manifestation of the law of large numbers 

q \ acting on the steady state probabilities. More importantly, because the model pinpoints mis- 

match repair as a basis of the rule, it is suitable for experimental verification/falsification. 
At this point, existing empirical survey supports the proposed mismatch-repair separation 
^ ■ between PR2 and non-PR2 genomes. 



O 



Introduction. Chargaff's parity rule with A = T, G = C for the nucleotide composition of DNA 
(flU) is mechanistically explained by Watson-Crick's double helix model of DNA ( |fT6ll ) because 
base A and T pair only with each other while base G and C pair only with each other along DNA's 
CN ■ complementary double strands. Chargaff and colleagues ([51 0) discovered a similar rule for single 

■ strands of DNA. For distinction, the former is referred to as Chargaff's Parity Rule (PR) and the 

latter as Chargaff's Second Parity Rule (PR2), respectively. Unlike PR, PR2 is not an exact rule 
but a statistical one. More specifically, according to a recent comprehensive test in [fT2|. PR2 holds 
for four of the five types of double stranded genomes: the archeal chromosomes, the bacterial 
^ . chromosomes, the eukaryotic chromosomes, most double stranded viral genomes, but it does not 

apply to the organellar (mitochondria and plastids) genomes nor to single stranded viral genomes 
or any type of RNA genome. Another study ([QQ|) also shows that the rule is length-dependent: for 
those that PR2 does apply, the longer the strand is the smaller the deviation from the rule becomes. 

In contrast to PR, the mechanistic basis of PR2 is still an open question (c.f. IU01 Base 
inversion, base substitution, base transposition are some of the suggested possible causes (lfTll8l[T5l 
IU). However, these suggested mechanisms fail to make clear and definitive distinction between 
those that obey the rule and those that do not. Various mathematical models are also proposed to 
explain the rule. However, they all led to a wrong prediction that single strands of a given length 
would eventually evolve to exhibit the exact rule after long enough time, which is inconsistent with 
the length-dependent property of the rule. This defect was noted in ( IfTOl ). 

Hence, any new model must make improvement in two fronts: its applicability only to the PR2 
genomes and its dependency on the genomic size. The mathematical model we propose here seems 
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to have accomplished both by limiting its applicability (a) to those organisms only with mismatch 
repair and (b) to the base evolution of any arbitrary single position along any single strand of the 
organism's DNA. The resulting model is a local Markov process for single bases whose steady 
state probabilities together with their manifestation via the law of large numbers give rise to the 
length-dependent empirical PR2. 

The Mathematical Model. Let m{B t — >• B r } denote the match or mismatch probability of a 
replicative base B r to a template base B t . The mathematical model of DNA replication with mis- 
match repair assumes the following. 

Mismatch Repair Error (MRE) Model: 

(a) DNA replication makes nucleotide mismatch errors. 

(b) At the moment of base replication, the mismatch error occurs independent of the 
nucleotide base position on the template single strand of DNA. 

(c) The match probabilities satisfy the following assumptions: 

i. The match probabilities between complementary bases are equal with re- 
spect to the complementary pairs: 

m{A — > T} = m{T — > A} and m{G — > C} = m{C — > G} 

ii. The match probability of a base to its complementary or itself is the same as 
that of the base's complementary to the complementary or itself: 

m{A — > A, or T} = m{T — > A, or T} and m{G — > G, or C} = m{C — > G, or C} 

iii. The conditional mismatch probability of a pyrimidine to its non-complementary 
purine equals that of its complementary purine to the purine's non-complementary 
pyrimidine: 

m{A -> C|A y4 T} = m{T -> G|T A} 

m{G -> T|G y4 C} = m{C -> A|C ~h G} 

(d) Mismatch repair makes strand recognition errors independent of complementary 
bases, and as a result it make base substitution errors. 

Hypothesis (a) is self-evident, which applies to all DNA, with or without mismatch repair. 
Hypothesis (b) is a localized and symmetrical assumption that each base position is as critical 
or ordinary as any other base position. It should be taken to be a primary approximation of this 
aspect of the replication. This means any assumption about global interactions along single strands 
of DNA, such as the stem-loop hypothesis of O, or codon position bias asymmetric assumption 
from [TTOl . can be taken as secondary approximations or corrections for future refinement of the 
model. 
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Hypotheses i,ii(c) implies that the mismatch probability of a base to itself is the same as that of 
the base's complementary to the complementary: 

m{A — > A} = m{T — > T} and m{G — > G} = m{C — > C} 

Alternatively, Hypothesis ii(c) can also be stated in terms of the mismatch conditional probabilities: 

m{A -> A|A y4 T} = m{T -> T|T y4 A} 

m{G -> G|G y4 C} = m{C -> C|C y4 G} 

In other words, Hypotheses i,ii(c) are made mainly on the complementarity of bases, which would 
be automatically true if DNA were a binary code in either the AT-pair or the GC-pair. 

On the other hand, however, Hypothesis iii(c) is made mainly on the steric characteristics of 
the pyrimidines and purines, i.e., when an A <A T (resp. G <A C) mismatch occurs, the conditional 
mismatch probability between A^C (resp. G^T) mismatch is the same as the T^G (resp. C— >A) 
mismatch. As a result of Hypotheses i,ii,iii(c), the remaining conditional mismatch probabilities 
are forced to satisfy: 

m{A — > G|A y4 T} = m{T — > C|T y4 A} 

m{G — > A|G y4 C} = m{C — > T|C G} () 

More details can be found from the match/mismatch probability diagrams of FigJTJ 

What separates our model from all others is Hypothesis (d). Although similar assumptions 
such as base inversion, inverted transposition, stem-loop substitution, codon position-bias were 
made for other models, the precise and systematic mechanisms that made such operations possible 
were either not known well enough or conjectured too broadly for all types of DNA, to which 
we know by the result of lfT2l that PR2 does not always apply. As a result, none is incorporated 
into the current model, nor is any repairing mechanism outside the phase of replication such as 
excision repair to spontaneous deamination of cytosine. In contrast, it appears that most types 
of double- stranded chromosomes are equipped with mismatch repair ( lfT3lO . More interestingly, 
mismatch repair in mitochondria of any organism is not known to exist according to lfl4ll . even 
though other types of repair mechanisms may exist (0) which do not negate the model under 
consideration. This distinction between DNA having or not having mismatch repair is consistent 
with the applicability of PR2 established in |[T2~1|. As a result of this hypothesis, our MRE model 
applies only to DNA types with mismatch repair satisfying these four hypotheses. 

The Result. From a modeling perspective, Hypotheses (c,d) imply that the MRE model is a 
Markov process model for any arbitrary position of any single strand of DNA since the transition 
probability at the position depends only on its current base. The Markov model is illustrated in 
FigfJJ referred to as the mismatch repair error chart, or MRE chart for short. Take, for example, 
a nucleotide A on a template single strand of DNA as shown at the top of the left chart. With a 
probability < 1 — a < 1, the process correctly replicates A's complementary base T, showing 
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__T~ -T- -A T-- --G T- -C- replicated base -C- --C-- --G C- -A C-- --T-- 

.-A- -A T- --A-- --C-- -A G- -G - --G-- --C-- -G T- --G-- --A-- 

Figure 1: MRE charts for base A and base G. Parameters a, d denote the mismatch probabilities 
m{A <A T},m{G <A C}, respectively, and bi,ei denote the conditional mismatch probabilities. 
For example, b 2 = m{A — > C|A y4 T} = m{T — > G|T A} and similarly, e 2 = m{G — > 



T|G t4 C} = m{C -> A|C t4 G}. 

down the left most branch, but incorrectly with < a < 1 probability. Of that fraction of mis- 
matches, for a fraction of < bi,b 2 each the process mismatches it with a G or a C, respectively, 
with bi + b 2 < 1. For the remaining \ — (bi + b 2 ) fraction of mismatches, an A is mismatched to 
the original A. These assumptions are represented by the middle three branches. 

For organisms which proof read and repair mismatches, the MRE chart continues one level 
down to the bottom branches. In such cases, Hypothesis (d) assumes that they do not always 
distinguish the replicative strand from the template strand 100% of the time, and make strand 
recognition errors independent of the replicating bases in question. Thus, for < c < 1 fraction 
of time, however insignificant or small it may be at this point, the process mistakes the replicative 
strand for the replicating template, and proceeds to turn the original base A into a base T if the 
mismatched pair is an A-A pair, with the hyphenated second base position being the template base. 
Similarly, with a probability < c < 1 each, a G-A pair is transformed to a G-C pair and a C-A 
pair to a C-G pair, all changing the template base A to a non-complementary base. With 1 — c 
probability, the process correctly identifies the template strand, and proceeds to keep the original 
base A in the template strand when the replication is completed. This particular illustration shows 
the case that the original mismatch error is corrected by a base T, which needs not to be the 
case. In fact, what is only assumed and important is that with the 1 — c probability the original 
base A is preserved along the template strand upon replication. The replicated new base along 
the complementary strand may or may not be the complementary base of the old base upon the 
completion of the mismatch repair. That is, the most general depiction of the chart would replace 
the base T by a place holder, O, for any of the four bases. 

Interchanging A and T and interchanging the parameters bi, b 2 in base A's MRE chart gives 
rise to base T's MRE chart if the replicating base at the top starts with a T. Because of Hypotheses 
(c,d), all the probability distributions for T's chart are the same as base A's chart except for the 
conditional mismatch probabilities of Hypothesis iii(c) and the consequence Eq.© of Hypothesis 
(c), which are accounted for by interchanging b\ and b 2 . The transition probabilities from base T 
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are listed on the second row of the transition probability matrix P below. 

Similarly, the diagram on the right is base G's MRE chart with probability distributions in 
d, ei, e 2 , f, which may not be respectively the same as a, bi, b 2 , c for A and T. In particular, a and 
d are very likely not to be the same, so are for bi and e^. According to Hypothesis (d), we have 
c = f. But there will be no changes to the result if we keep c, / distinct. Also, interchanging G 
and C in G's MRE chart gives rise to C's MRE chart. 

The mismatch repair charts allow us to derive the transition probabilities between bases one 
replicating generation a time. More specifically, take the A's MRE chart for an example. Let 
p{A — > A} be the probability that a base A from a template strand remains an A upon replication. 
Similarly, let p{A — > T} be the the probability that a base A from the template strand is substituted 
by a T upon replication, and similar notation applies to p{A — > G},p{A — > C}. Then these 
probabilities can be tabulated from probability distributions from the MRE chart as follows: 

p{A -> T} = a(l - (61 + b 2 ))c 
p{A —> G} = ab\C 
p{A — > C} = ab 2 c 

p{A -»• A} = 1 - p{A -> T} - p{A -> G} - p{A -> C} = 1 - ac 

The first expression, for example, is obtained by following from the top of the chart down the direct 
branches leading to the substituting base T and multiplying all the probabilities along the branches. 
The same for the transition probabilities from A to G and C respectively. The probability from 
A to A can either be obtained by the formula above or by summing all probabilities of the A-to- 
A branches (four branches in all) of the chart. Exactly the same tabulation gives the transition 
probabilities for all other bases. 

Using matrix entry notation, we denote 

Pu = p{A -> A}, Pl2 = p{A -> T}, Pl3 = p{A G}, p 14 = p{A C}. 

That is, the row and column numbers, 1, 2, 3, 4, are in correspondence with A, T, G, C, respec- 
tively. Similar notation extends to other bases as well. As a result we obtain the following transition 
probability matrix for our MRE model: 



1 — ac a(l — (61 + b 2 ))c ab^c ab 2 c 

a(l — {b\ + b 2 ))c 1 — ac ab 2 c ab\c 

deif de 2 f 1 - df d(l - (e 1 + e 2 ))f 

de 2 f de 1 f d(l - (e 1 + e 2 ))f 1 - df 



What follows are standard textbook properties of transition probability matrixes of Markov chains 

(c.f.ini). 

1 . Denote the nth iterate of the transition matrix by 

P n = [pg°]- 
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(n) 

Then we know that it is again a transition probability matrix with its entry p\- ' represent- 
ing the probability of a base % becoming a base j after the nth generation of replication. 

2. Because the transition matrix P has all positive entries (hence is irreducible and er- 
godic), by the Perron-Frobenius Theorem, the limit lim n ^oo p\j exists and the limit is 
independent of the initial base i: Hm^oopjj = 7Tj, satisfying < ttj < 1 and 

71~1 + 7T 2 + 7T 3 + 7T 4 = 1. 

In terms of DNA replication with mismatch repair, probability ttj is the steady state 
probability of finding base j at any base position of any single strand of DNA. This 
means, regardless of the initial base i at that position, after sufficiently many generations 
of replication, the probability of finding base j at the position is iij. We denote the four 
steady state probabilities by 7Ti = p{A}, 7r 2 = p{T}, tt 3 = p{G}, ir 4 = p{C}. 

3. The transition matrix P has Ai = 1 to be the largest eigenvalue in magnitude and it is 
simple. Moreover, the steady state probability vector n = [ttx, ir 2 , 7r 3 , 7r 4 ] is the only left 
eigenvector not counting scalar multiple (or the transpose ix 1 is the only right eigenvector 
of eigenvalue 1 for the transpose matrix P*). Because of this property, we can use the 
eigenvector equation irP = n to find ir explicitly as 

= . 77 • t~t • ~r, ■ - 7 7[d(e 1 +e 2 )f\ d(e l +e 2 )f, a{b l +b 2 )c, 0(61+62)0], 

2(a(6i + b 2 )c + d(ei + e 2 )f) 

with equal probabilities 7Ti = ir 2 , tt 3 = 7r 4 for complementary bases. This, to be ex- 
plained later, is the basis for the empirical law of PR2. 

4. By the Perron-Frobenius Theorem, the remaining eigenvalues of P are less than 1 in 
magnitude. In this case, they can be explicitly found as: 

A 2 = 1 - a(b x + b 2 )c - d(e 1 + e 2 )f 

A 3i4 = 1 - ac — df + |(a(6i + b 2 )c + d(e x + e 2 )f) ± IVD, where 
D = (-2ac + 2df + a(6i + b 2 )c - d{e 1 + e 2 )ff 
+4a(6i - 6 2 )cd(ei - e 2 )f 

It has the property that for any initial probability distribution q = [qi,q 2 ,q 3 , g 4 ] with 
<7i + <?2 + <?3 + 94 = 1, qP n converges to n at a rate no greater than the order of \fi\ n with 
fj, = max 2 <j<4{|Aj|} < 1. This gives a temporal estimate for the convergence rate to the 
steady state probability at the nth generation of replication. 

We now conclude that Hypothesis (b), Property 2, and the law of large numbers imply the empirical 
PR2. More specifically, let Fi{i} be the length-averaged frequency of base i, called the ensemble 
frequency, then we must have lim^oo -FlW = p{i} for i = A, T, G, C. To see this, we borrow 
a prototypical explanation from coin tossing: The steady state probabilities of the head and tail 
of a coin can be approximated by tossing an ensemble of L many identical coins and counting 
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Figure 2: Homo sapiens chromosome 22 genomic contig ref|NT_l 13818.1 |Hs22_l 1 1537: 1- 
17927 (from Human Genome Project) with length L = 17,927 and ensemble frequency Fl = 
(0.1993,0.2007,0.2917,0.3083). (a) In the 5' -> 3' direction, the sequence of ensemble fre- 
quencies F N with length N = 1, 2, . . . , L is generated. The convergence error sequence, 
X\=a T G c \Fn(i) — F L (i)\, is plotted against N together with its best fit to the curve a + b/ y/~N. 
(b) The same plot except for the opposite 3' — > 5' direction, (c) The same plot except (1) the data is 
averaged over 10 runs and (2) each run is done for a random permutation of the original sequence. 



the ensemble frequencies in which case the greater the assemble size L is the closer the ensemble 
frequency to the steady state probability distribution becomes. In fact, according to the law of 
large numbers (or the Central Limit Theorem), the convergence rate has the order of 1 jyL. Figf2] 
gives a typical simulation of the convergence for a short contig of the Human chromosome 22. 
Simulations with longer sequences are not shown here because in general the longer a sequence 
the better the convergence fit becomes. 

Discussion. It is interesting to note a contrast between mismatch repair and other possible or con- 
jectured processes such as base inversion, transposition, stem-loop substitution, codon position- 
bias asymmetry. For organisms without mismatch repair, the MRE chart stops at the first or the 
second branches with the template base well-entrenched in its position. It is not clear what mech- 
anistic operations these processes adopt to systematically substitute the template base. If such 
processes only subject the template base to a sporadic, disorganized chance substitution as implied 
when such systematic mechanisms are lacking, it is hard to see how a regular and steady- state-like 
pattern such as PR2 can arise. It is also important to note that the model starts with a single strand 
separated from a double stranded mother DNA during replication and ends with the completion of 
replication. This means that any damage or repair that may take place before or after replication 
is not relevant to the model at its current formulation. In other words, if such types of damages 
and repairs are not the results of some systematic, repetitive, and normal cellular functions, but 
rather disorganized random events, then the lack of regular, steady- state-like patterns in non-PR2 



7 



genome is not a contradiction to the model. 

On the surface of it, mismatch repair enhances replication fidelity by enforcing PR. However, 
if the proposed model is correct, it will imply that mismatch repair results in PR2. The question 
is then: is PR2 also some sort of genomic fidelity that needs to be preserved? If so, what is the 
evolutionary advantage of having it? A probable answer can be found in the communication theory 
for DNA replication proposed in Q and it will be addressed elsewhere. 

Past modeling attempts were not successful mathematically because they did not make the 
critical distinction between the steady state probability at individual base position and the ensemble 
frequency of whole length along single strands of DNA. To a lesser extent, they failed because they 
did not limit their scopes to different types of DNA by pinpointing the mechanistic processes which 
eventually but systematically led to the manifestation of PR2 we observe today. These assessments 
motivated the model proposed here which should be viewed as a first attempt to address both issues 
simultaneously. Future refinements of the model should be done by modifying and/or expanding 
the MRE hypotheses (a,b,c,d) in a manner as concise as we did here. 

Non-PR2 genomes can be used as an empirical control group to guard against possible placebo 
effect of any proposed mechanism for PR2. For example, the hypothesis of base inversion (|fT|) 
can be ruled out to be the cause of PR2 because it applies to both PR2 and non-PR2 genomes. 
Our proposed mismatch-repair mechanism, on the other hand, has so far withstood this scrutiny by 
keeping the empirical separation between PR2 and non-PR2 genomes. 

Because of its conciseness, the model is suitable for experimental test. More specifically, an 
experiment can start with many identical double strands of DNA about 6,000 or more bp long 
(c.f. Fig© with an ensemble frequency disobeying PR2. Divide them into two groups. Repli- 
cate one group by a replicator with mismatch repair and replicate the other group by a replicator 
without mismatch repair. After a sufficiently long period of time, calculate the averaged ensemble 
frequency from each group and compare. The experiment would be supportive of the model if 
the ensemble frequency of the mismatch repairing replicator is closer to PR2 than that of the null 
replicator. The model can also be rejected if one can find a non-PR2 genome with mismatch repair. 
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